function [ cftimage ] = stoltmig(cftdata, fs, dx, ERMv, gamma)
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here
[nt, nx] = size(cftdata);
f = (-nt/2:nt/2-1)*fs/nt;
kx = (-nx/2:nx/2-1)/dx/nx;
[kx,f] = meshgrid(kx,f);
fkz = ERMv*sign(f).*(sqrt(kx.^2 + f.^2/ERMv^2)*cos(gamma) + kx*sin(gamma));
cftimage = interp2(kx, f, cftdata, kx, fkz, 'linear', 0);
end

